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Abstract 

The physics of the tv phase shift in ferromagnetic Josephson junctions 
jrt ■ enables a range of applications for spin-electronic devices and quantum 

computing. In this respect our research is devoted to the evaluation of 
the best materials for the development and the realization of the quan- 
'"^ ■ turn devices based on superconductors and at the same point towards the 

Ch ' reduction of the size of the employed heterostructures towards and below 

O , nano-scale. In this chapter we report our investigation of transitions from 

, >^/ to TT states in Nb Josephson junctions with strongly ferromagnetic bar- 

riers of Co, Ni, Ni8oFe2o (Py) and Fe. We show that it is possible to fabri- 
cate nanostructured Nb/ Ni(Co, Py, Fe)/Nb 7r-junctions with a nano-scale 
^ ' magnetic dead layer and with a high level of control over the ferromag- 

^\ I netic barrier thickness variation. In agreement with the theoretical model 

^T ■ we estimate, from the oscillations of the critical current as function of the 

00 ' ferromagnetic barrier thickness, the exchange energy of the ferromagnetic 

>— -^ , material and we obtain that it is close to bulk ferromagnetic materials 

[■■*-. ■ implying that the ferromagnet is clean and S/F roughness is minimal. 

c7^ ' We conclude that S/F/S Josephson junctions are viable structures in the 

00 , development of superconductor-based quantum electronic devices; in par- 

^^ ' ticular Nb/Co/Nb and Nb/Fe/Nb multilayers with their low value of the 

magnetic dead layer and high value of the exchange energy can readily 
be used in controllable two- level quantum information systems. In this 
^\^ ' respect, we discuss applications of our nano-junctions to engineering mag- 

JH , netoresistive devices such as programmable pseudo-spin-valve Josephson 

structures. 



1 Introduction 

Conventional electronics is based on the transport of electrical charge carriers, 
however the necessity to have more versatile and efficient devices on micro- and 
nano-scale has diverted the attention towards the spin of the electron rather than 
its charge: this is the essence of spintronics. Although spintronic devices were 
traditionally built with semiconductors and/or ferromagnetic materials, nowa- 
days hybrid structures constituted by superconductors (S) and ferromagnetic 
(F) materials are emerging as promising alternatives. They serve as candidates 
to realize and investigate systems useful both as testgrounds to probe fundamen- 
tal physics questions, and as realistic elements for application in the spintronics 
industry thanks to the improved degree of control on the spin polarized current. 

In hybrid S/F structures it has been verified that, due to the simultaneous 
presence of these two competitive orders, peculiar effects appear: the supercon- 
ductivity is reduced by the spin polarization of the F layer, while due to the 
proximity effect the Cooper pairs can enter in the superconductor resulting in a 
oscillatory behavior of the density of states. The last effect opens the doors to- 
wards new and exciting developments for the realization of quantum electronic 
devices. In fact in S/F/S systems, due to the oscillations of the order parameter, 
some oscillations manifest in the Josephson critical current as a function of the 
ferromagnetic barrier thickness evidencing the presence of two different states, 
and TT, corresponding to the sign change of the Josephson critical current. 

Our research, presented in this chapter, fits in this rapidly developing and 
exciting area of condensed matter physics [T]. Within the context of the spin 
polarized devices based on the interplay between superconductivity and fer- 
romagnetism, we have fabricated S/F/S nano-structured Josephson junctions 
constituted by a low temperature superconductor. Niobium (Nb), and strong 
ferromagnetic metal, Nickel(Ni), NigoFe2o (Py), Cobalt (Co) and Iron (Fe), and 
we have investigated their magnetic and electrical properties. These structures 
have evidenced a small magnetic dead layer and oscillations of the critical cur- 
rent as a function of the ferromagnetic barrier. These oscillations show excellent 
fits to existing theoretical models. We also determine the Curie temperature 
for Ni, Py and Co. In the case of Co and Fe we estimate the mean free path 
to confirm that the oscillations are in the clean limit and from the temperature 
dependence of the IqRn product we show that its decay rate exhibits a non- 
monotonic oscillatory behavior with the ferromagnetic barrier thickness. We 
investigate the presence of the Shapiro steps on the I vs V curve applying a 
microwave and the effect of the magnetic field on the maximum supercurrent. 
In this last case we find oscillations of the maximum of the supercurrent cor- 
risponding to a Fraunhofer pattern. Finally, in the case of Co barrier, focusing 
in detail on a single O-tt phase transition we show evidence for the appearance of 
a second harmonic in the current-phase relation at the minimum of the critical 
current. For the high value of the exchange energies and small magnetic dead 
layer the S/F/S structures with Co or Fe barrier can be considered as good 
candidates for the realization of quantum devices. 
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Figure 1: Andreev reflection process for a metal with spin polarization P = 
(a) and P = 100% (c); experimental measurements of a N/S interface (b) and 
F/S junction (d). Figure adapted from Ref.[3]. 



2 F/S Junctions: Basic Aspects 

The aim of this section is to give an outline of the physics of the Ferromagnetic- 
Superconducting (F/S) interfaces (for a review see [2])- Andreev reflection plays 
an important role to understand the transport process in F/S junctions. The 
Andreev reflection near the Fermi level preserves energy and momentum but 
does not preserve spin, in other words the incoming electron and the reflected 
hole have opposite spin. This is irrelevant for the transport in N/S junctions 
(N stands for normal metal) due to the spin-rotation symmetry. On the other 
hand at F/S interfaces, since the spin-up and spin-down bands in F are different, 
the spin flipping changes the conductance profile. In particular in fully spin- 
polarized metals all carriers have the same spin and the Andreev reflection is 
totaly suppressed; therefore, at zero bias voltage the normalized conductance 
becomes zero (see Fig. [1]). In general, for arbitrary polarization P, it can easily 
be shown that 

^=2(1-P) (1) 

where, in terms of the spin-up N-^ and spin-down N^ electrons, the spin polar- 
ization is defined as 

{Nt + n^Y ^^ 

From this analysis it has been shown that the point contact measurements 
based on the Andreev reflection process give a quantitative estimation of the 




Figure 2: Oscillatory behavior of the exponential decay of the superconducting 
order parameter at the F/S interface. 



polarization of the ferromagnetic material [21 E] ■ In fact from the reduction of 
the zero bias conductance peak it is possible with a modified Blonder-Tinkham- 
Klapwijk model to estimate the polarization P. 

In a ferromagnet in proximity to a superconductor, the coherence length S,f, 
due to the presence of the exchange field E^x , is given by: 
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2{Tr KbT + lEex) ' 



(3) 



where D is the diffusive coefficient. If the exchange energy is large compared to 
the temperature, Eex > ksT, the coherence length is much shorter than in the 
case of the N/S proximity effect. In addition to the reduced coherence length 
compared to typical N/S structures, a second characteristic property arises from 
the complex nature of the coherence length: the induced pair amplitude oscil- 
lates spatially in the ferromagnetic metal as a consequence of the exchange field 
acting upon the spins of the two electrons forming a Cooper-pair [5l [6l |7l [8] 
(see Fig. [5]). This oscillation includes a change of sign and by using appropriate 
values for the exchange energy and F-layer thickness, negative coupling can be 
realized. We refer to the state corresponding to a positive sign of the real part 
of the order parameter as "0-state" and that corresponding to a negative sign 
of the order parameter as "7r-state" . 

Demler et al. [9] gave a qualitative picture of this oscillatory character. They 
considered a Cooper pair transported adiabatically across an F/S interface, with 
its electron momenta aligned with the interface normal direction. The pair 
entering in the F region decays exponentially on the length scale of the normal 
metal coherence length. Then the up-spin electron, oriented along the exchange 
field, decreases its energy by ft, = Ef,x/H, where E^x is the exchange energy of the 
F layer. On the other hand, the down-spin electron increases its energy by E^^x- 
To compensate this energy variation, the up-spin electron increases its kinetic 
energy, while the down-spin electron decreases its kinetic energy. As a result the 
Cooper pair acquires a ccnter-of-mass momentum Q — 2Ef,x/vFi which implies 
the modulation of the order parameter with period t^vf / Ef.x, where vp is the 
Fermi velocity. 




Figure 3: Schematic of the Andreev reflection process in a F/S junction, the 
momentum shift, App, is dominated by the spin sphtting of the up and down 
bands. 



As a consequence of the oscillations of the order parameter, a similar os- 
cillatory behavior is observed for the density of states. This behavior can be 
explained by considering the spin effect on the mechanism of Andreev reflections 

The process is illustrated in Fig. [3] using the energy-momentum dispersion 
law of the normal metal: in the case of a N/S interface an incoming electron 
in a normal metal N with energy lower than the superconducting energy gap A 
from the Fermi level can be reflected into a hole at the N/S interface (Andreev 
Reflection) [TUj. If the normal layer is very thin the density of states in N is 
close to that of the Cooper pair reservoir. The situation is strongly modified if 
the normal metal is ferromagnetic. As Andreev reflections invert spin-up into 
spin-down quasiparticles and vice versa, the total momentum difference includes 
the spin splitting of the conduction band: APp ~ Q. The density of states is 
modified in a thin layer on the order of ^p. In particular, the interference 
between the electron and hole wave functions produces an oscillating term in 
the superconducting density of states with period Eex/fivp. This effect has been 
observed experimentally by Kontos et al. [11] 



2.1 Theory of the Josephson yr-junctions 

Due to the spatially oscillating induced pair amplitude in F/S proximity struc- 
tures it is possible to realize negative coupling of two superconductors across 
a ferromagnetic weak link (S/F/S Josephson junctions). In this case of nega- 
tive coupling, the critical current across the junction is reversed when compared 
to the normal case giving rise to an inverted current-phase relation. Because 
they are characterized by an intrinsic phase shift of tt, these junctions are called 
Josephson it -junctions. 

One of the manifestations of the tt phase is a non-monotonic variation of 
the critical temperature and the critical current, with the variation of the fer- 
romagnetic layer thickness {tp) [2]. Referring to the current-phase relation 



Ic = Iq sin (j>, for a S/I/S Josephson junction the constant /q > and the mini- 
mum energy is obtained for (j> = 0. In the case of a S/F/S Josephson junctions 
the constant Iq can change its sign from positive to negative indicating the 
transition from the 0-state to 7r-state. Physically the changing in sign of /q is 
a consequence of a phase change in the electron pair wave function induced in 
the F layer by the proximity effect. Experimentally, measurements of Ic are in- 
sensitive to the sign of /q hence the absolute value of Iq is measured, so we can 
reveal a non-monotonic behavior of the critical current as a function of the F 
layer thickness. The vanishing of critical current marks the transition from to 
TT state. The dependence of the critical current on the thickness of the ferromag- 
netic layer in S/I/F/S junctions has been experimentally investigated by Kontos 
et al. [12]. The quantitative analysis of the S/F/S junctions is rather compli- 
cated, because the ferromagnetic layer can modify the superconductivity at the 
F/S interface. Then other parameters, such as the boundary transparency, the 
electron mean free path, the magnetic scattering, etc can affect the critical cur- 
rent. It is outside the purpose of this chapter to derive explicitly the expression 
of the critical current as a function of the F layer, for the interested reader we 
refer to this review [2]. 

The majority of experimental studies have concentrated on weak ferromag- 
nets where Ef,^ ~ KbTc, where T^ is the superconducting critical temperature, 
resulting in multiple oscillations in /j, with temperature and tp- In the case of 
strong ferromagnets, where Tcurio ^ Tcj only oscillations of Ic with tp, and not 
with temperature, are observed. In this chapter we will present the study of 
the oscillations of the critical current as a function of tp for S/F/S Josephson 
junctions with strong ferromagnetic barriers. 

The generic expression of the critical current as a function of F layer is given 
by: 
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where ii is the thickness of the ferromagnet corresponding to the first minimum 
and IcRN{to) is the first experimental value of IcRn {Rn is the normal state 
resistance), and fi and ^2 are the two fitting parameters. Eq. granges in the 
clean and in the dirty limit. In particular, in clean limit, tp < L where L 
is the mean free path of the F layer, ^2 = vph/2Ecx- In this way, known ^2 
and estimating the Fermi velocity from reported values in literature, one can 
calculate the exchange energy of the ferromagnetic barrier. 

Then, in the case of clean limit the oscillations of IcRn vs tp can be modeled 
by a simpler theoretical model |T?] given by: 

IcRNitp) - IcRNito) 2Ecxt^/hvf ' ^^> 

where in this case the two fitting parameters are vp and Eg^. 

On the other hand, in dirty limit tp > L we can model the oscillations [14] 



by the following formula: 
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where A is the superconducting order parameter, LOm is the Matsubara frequency 
and is given by Um = '!TTkB{'2m + 1) where T is the transmission coefficient and 
m is an integer number, k^^ = {1 + 2 \ w^ | t/?*) — ^lE^.^r/h and fj, — cos where 
6 is the angle the momentum vector makes relative to the distance normal to 
the F/S interface. L is given by vpT and t is the momentum relaxation time. 
In this case the fitting parameters are vp, E^x and the mean free path L of the 
ferromagnetic layer. 

3 Josephson Junction Fabrication 
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Figure 4: A schematic picture of the S/F/S heterostructure used in this chapter. 



To investigate the physics of the ferromagnetic 7r-junctions and their possible 
applications for spin-electronic devices we have realized and characterized Nb 
(250 nm thick) / Ferromagnetic layer / Nb (250 nm thick) Josephson junctions. 
As ferromagnetic layer we have used: Ni, Py, Co and Fe. In the case of Co, Fe 
and Py their respective thicknesses t-p were varied from ~ 0.5-5.5 nm while, in 
the case of Ni, ^ni was varied from ~ 1.0 to 10 nm. To assist processing in a 
focused ion beam (FIB) microscope, a 20 nm normal metal interlayer of Copper 
(Cu) was deposited inside the outer Nb electrodes, but 50 nm away from the 
Fe barrier. We remark that the 20 nm of Cu is a thickness smaller than its 



coherence length, so it is completely proximitized into the Nb, and it does not 
affect the transport properties of the Josephson junction. Refer to Fig. [4] for an 
illustration of our he □ 
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Figure 5: An illustration of the sputtering process: the ferromagnetic interlayer 
thickness is varied as a function of the angle 9 from the pre-sputter position. 

The heterostructures employed in this chapter, have been deposited by d.c. 
magnetron sputtering. In a single deposition run, multiple silicon substrates 
were placed on a rotating holder which passed in turn under three magnetrons: 
Nb, Cu and the ferromagnetic layer (Fe, Co, Py or Ni). The speed of rotation 
was controlled by a computer operated stepper motor with a precision angle of 
better than 3.6° and each sample was separated by an angle of at least 10°. Prior 
to the deposition each target material was calibrated and the deposition rates 
were estimated with an Atomic Force Microscopy (see table [1] for a summary 
of the deposition parameters for all target materials presented in this chapter). 
In the case of the ferromagnetic materials (Co, Py, Fe and Ni) the rates of 
deposition, and hence ip, were obtained by varying the speed of each single chip 
which moves under the ferromagnetic target while maintaining constant power 
to the magnetron targets and Ar pressure. This was achieved by knowing the 
chip position relative to the target material (9) and by programming the rotating 
flange such that a linear variation of ferromagnetic thickness with 0, d{tp)/d9, 
was achieved, tp is inversely proportional to the speed of deposition, tp oc l/Vj, 
and hence it can be shown that in order to achieve a linear variation of tp{6) 
one programmes the instantaneous speed, at position 9 and time t seconds (i.e. 
V{9)t), of deposition according to 

where Vi is the initial speed and Vf is the final speed in units of rpm. This 
method of varying tp guaranteed, in all cases, that the interfaces between each 
layer were prepared under the same conditions while providing precise control of 
the F layers [TS]. Fig. [5] shows an illustration of the sputtering process. For each 



Table 1: A summary of the deposition parameters for all materials 
sputtered, tp refers to the expected film thickness. 



Target material 



Rate 



Power Speed range 



{nm/W at 1 rpm) (W) 



(rpm) 



tp 
(±0.2 nm) 



Nb 
Cu 
Ni 
Co 
Py 

Fe 



6.89 


90 


.028 


250 


4.39 


30 


0.22 


20 


2.12 


40 


0.20-4.2 


0.5- 10.5 


1.8 


40 


0.36-3.6 


0.5 - 5.0 


1.64 


40 


0.33-3.3 


0.5 - 5.0 


2.0 


40 


0.22-2.2 


0.5-5.0 



run, simultaneously, our heterostructures have been deposited on 5 x 5 mm^ and 
5 X 10 mm^ Si02 substrates. The first ones have been used for the magnetic 
measurements, the second ones for the realization of Josephson junctions. 

3.1 X-ray Measurements 



o Expected 

• Observed (X-ray) 
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Figure 6: tco(ohserDed) plotted with tco(expected) extracted from simulation data, 
in the inset: low angle x-ray data plotted with the equivalent simulation for 
sample 1 where tco{expected) =5 nm and tco(observed) =5±0.2 nm. The mean 
deviation is estimated to be ±0.2 nm. 



To confirm our precise control over the ferromagnetic thickness variation we 
performed low angle X-ray reflectivity of a set of calibration Nb/Co/Nb thin 
films where the Nb layers had a thickness of 5 nm and the Co barrier thickness 



was varied from 0.5 nm to 5.0 nm [TS]. A series of low angle X-ray scans were 
made and the thickness of the Co layer (ico(o6ser«ed)) was extracted by fitting 
the period of the Kiessig fringes using a simulation package. It was found 
that our expected Co barrier thicknesses, ico(e£cpected)j was well correlated with 
tco(observed) with a mean deviation of ±0.2 nm. In the inset of Fig. [S]we show 
an example of the low angle x-ray data plotted with the equivalent simulation 
data where tco(obseri;ed) was extracted, while in Fig. [6] we report a comparison 
of tcoiobserved] with tco{expected) and its deviation. 



4 Nanoscale Device Process 



Contact pads 




Figure 7: (a) "Cam 39" mask design for a 10 x 5 mm^ chip used in this chapter, 
showing contact pads for wirebonding. (b) Central detail: on each chip we 
realize 14 possible junctions with thinnest track widths of 4/im (c). 

We can summarize the realization of the Josephson junctions in three dif- 
ferent steps: (i) patterning the films using optical lithography to define the 
micron scale tracks and the contact pads. Our mask permits etching of at least 
14 devices, which allows us to measure numerous devices and to derive good 
estimates of important parameters, like, for example, characteristic voltage (see 



10 



Fig. [7| to see an illustration of the mask); (ii) broad beam Ar ion milling (3 
mAcm"^, 500 V beam) to remove unwanted material from around the mask 
pattern, thus leaving 4/im tracks for subsequent FIB work; (iii) FIB etching of 
micron scale tracks to achieve vertical transport [TH [T71 [TS] . 

In particular to realize our devices we have used a three-dimensional tech- 
nique [IH]- The wedge holder used in this chapter was designed by D.-J. Kang 
and it is schematically shown in Fig. [H it is constituted by three sample lodg- 
ings, one in horizontal and two at 45°. Once loaded the sample on the 45° 
lodging, we can rotate the stage of the FIB at 45°, in this way the beam is 
perpendicular to the surface of the sample, and the first cut is done; then the 
sample holder is rotated of 180° around an axis normal to the sample stage, to 
permit the vertical cut to be done (see Fig. [8]). This setup allows to load two 
samples for each run. 

(a) 



sample 





Figure 8: Schematic of the sample holder used in FIB (a) showing two axes of 
rotation, 45°, to horizontal milling (b) and 180° to vertical milling. 

The fabrication procedure is shown in Fig. [9] [19]. A first box with area 
4 X 2fim^ is milled with 150 pA to realize tracks of about 700 nm (b). The time 
of milling is about 1 — 3 minutes, and this milling can be calibrated using the 
stage-current /end-point detection measurement. Fig. [10] shows how the milling 
of different layers can be distinguished, the first peak corresponds to the Nb, 
the two small peaks correspond to the two Cu layers, and finally the intensity 
decreases approaching to the Si02 substrate. The sidewalls of the narrowed 
track are then cleaned with a beam current of 11 pA. This removes excessive 
gallium implantation from the larger beam size of the higher beam currents. 
The cleaning takes ^ 25-30 seconds per device. The track width is then < 500 
nm. The sample is then tilted to 9 = 90°, and the two cuts are made with a 
beam current of 11 pA to give the final device (Fig. [9] (d)) with a device area in 
the range of 0.2 — 1 /^m^. This technique permits to achieve vertical transport 
of the current: in Fig. [5] (e) we show a schematic picture of the current path. 
In Fig. [Tl]thc final FIB image of a Nb/Cu/Co/Cu/Nb device is presented. 
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Figure 9: FIB procedure for device fabrication: the initial trilayer after the 
photolithography and ion-milling (a) is cut with 150 pA beam (b) and then 
with 11 pA beam (c-d). (e) At the end two side cuts are realized to create the 
final device structure with a device area in the range of 0.2 — 1 /xm^ achieving 
vertical transport of the current. (Figure adapted from ref. [10]). 

5 Magnetic Measurements 

In this section we report magnetic measurements of Nb Josephson junctions 
with strongly ferromagnetic (F) barriers: Ni, Ni8oFe2o (Py), Co and Fe. From 
measurements of the magnetization saturation (Ms) as a function of the F thick- 
ness, our heterostructures have shown a magnetic dead layer ranging between 
0.5 nm and 1.7 nm. Then we give an estimation of the Curie temperature of 
the ferromagnetic layer. 

5.1 Measurement of the Magnetic Dead Layer 

In this section we explain the magnetic properties of Nb-Ni-Nb, Nb-Py-Nb, Nb- 
Co-Nb and Nb-Fe-Nb heterostructures as a function of F thickness [THl [HI [HI 

[m. 

To this aim wc have studied, using a VSM at room temperature, the hys- 
teresis loop of our heterostructures in order to follow the evolution of the mag- 
netization as a function of the applied magnetic field. In Fig. [T2| (inset). Fig. 
fT3l (a), [HI (a) and (b), we show a collection of hysteresis loops for different 
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Figure 10: End Point Detector. The graph shows the stage current as a func- 
tion of the miUing time, we can distingue the Nb layer, the two Cu layers, 
the Ferromagnetic barrier and when the current decreases we have reached the 
substrate. 




Figure 11: FIB image of a Nb/Cu/Co/Cu/Nb device, the two light gray lines 
are the Cu markers. 



thicknesses of the F barrier. We notice that both the Ms and the width of the 
hysteresis loop, as expected, decrease with decreasing F layer thickness, and the 
ferromagnetic order disappears when the F barrier goes to zero. Furthermore 
from the hysteresis loops we have measured the saturation magnetization as a 
function of the F barrier to extrapolate the magnetic dead layer. 
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• Experimental data 
Least-squares fit 
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Figure 12: Saturation magnetization versus thickness of the Ni barrier, inset: 
hysteresis loop for iwi = 4 nm. 
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Figure 13: (a) Hysteresis loops for different Py thicknesses, (b) Saturation 
magnetization versus thickness of the Py barrier. 



The presence of a magnetic dead layer has been reported in other studies of 
S/F/S heterostructures [22l[23l[15] and it can be explained as a loss of magnetic 
moment of the heterostructures. The magnetic dead layer can be due to numer- 
ous factors, as, for example, lattice mismatch at the Nb-F interface resulting in 
a reduction in the ferromagnetic atomic volume [24] and crystal structure which 
leads to a reduction in the exchange interaction between neighboring atoms. 
This loss in exchange interaction manifests itself as a loss in magnetic moment 
due to a collapse in the regular arrangement of electron spin and magnetic mo- 
ment and can lead to the suppression in Tcurio and Eex [lH [26l [27l [28] . Another 
factor can be the inter-diffusion of the ferromagnetic atoms into the Nb. Like 
in the case of lattice mismatch this would result in a breakdown of the crystal 
structure at the interface leading to a reduction in the exchange interaction. 
The knowledge of the magnetic dead layer is a crucial point within the imple- 
mentation of TT-technology, to guarantee the reproducibility and the control of 
the devices. Extrapolating with a linear fit the saturation magnetization as a 
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Figure 14: Hysteresis loops for different Co (a) and Fe (b) thicknesses. Satura- 
tion magnetization versus thickness of the Co (c) and Fe (d) barrier. From the 
linear fit (gray line) the magnetic dead layer is extrapolated. 



function of F thickness, we have obtained the estimation of the magnetic dead 
layer: ^1.7 nm for Ni (see Fig. [T^ . ^ 0.5 nm for Py (see Fig. [12] (b)) and 
- 0.75 nm for Co (see Fig. [H (c)). 

In the case of Fe the saturation magnetization was measured for three dif- 
ferent deposition runs (see Fig. [14] (d)). For each deposition we have obtained 
similar saturation magnetization and by extrapolating the least-squares fit of 
these data we have estimated a Fe magnetic dead layer of ~ 1.1 nm |T7]. In 
Fig. [T5]we summarize the magnetic moment vs thickness for Ni, Py, Co and Fe. 
From these data we can extrapolate the slope of Ms vs tp. From the theoretical 
model [29j_30^ the predicted slopes are: ~ 0.60 emu/cm^ for Ni, ~ 0.52 emu/cm'^ 
for Py, ~ 1.42 emu/cm^ for Co and ~ 2.6 emu/cm'^ for Fe. On the other hand 
from our experimental data the slopes of Mg vs ip for these ferromagnetic mate- 
rials are suppressed from these expected bulk values: 0.27 emu/cm'^ for Ni, 0.27 
emu/cm^ for Py, 1.0 emu/cm'^ for Co and 0.83 emu/cm'^ for Fe. We can argue 
that, in our case, we are not considering bulk materials, as reported from Slater 
and PauHng [29l[30j, instead our systems are constituted by F/S sandwiches, 
so the presence of a superconducting layer, Nb, can induce a weakening of the 
ferromagnetic properties of the F layer and a possible diffusion of the Nb into 
the F barrier. 



15 



6n 






A 


Co 


5- 




/ 


• 


Py 
Ni 


4- 




«y« /» 


* 


Fe 


3- 




/ y 






2- 


/ 


/ 




^^^^ 




/a 


7* 






1- 


/y^ 


r*^^^^^^^^ 








JJ><^ 


^^""^ 






0- 


^•^^^ 





1 ' 1 ' 


, . , 



1 



10 11 



tp (nm) 
Figure 15: Saturation magnetization vs Ni, Py, Co and Fe thickness at T= 300K. 

5.2 Calculation of the Curie temperature 

We have measured the thermal variation olthe saturation magnetization, M{T), 
of Co, Ni, and Py when sandwiched between thick Nb layers [15]. To be certain, 
M{T) with temperature was not weakened by a thermally activated diffusion of 
ferromagnetic atoms into Nb, or vice versa; at the interface we measured both 
the magnetization when warming and cooling. The warming and cooling data 
agreed for all three barrier systems up to a temperature of 620 K. Above this 
temperature, Ms was found to drop by virtue of thermally activated diffusion. 
We have modeled the warming and cooling data of Ms (lOK < T < 620K) with 
the following formula 

M(T)/M(0) = (1 - T/Tcurief, 

where M(0) is the saturation magnetization at absolute K, T is the measuring 
temperature, and /3 and Tcurio are fitting parameters. This gives Tcurie values 
of 1200 K for Co, 571 K for Ni, and 800 K for Py (see FigHH) in agreement 
with the bulk values. Data for Ni are the most reliable because we have a full 
data set. However, in any case these measurements provide that, in our metallic 
systems, interdiffusion at the ferromagnetic surface cannot be ruled out because 
the ferromagnetic layer is known to form a variety of magnetic and nonmagnetic 
alloys with Nb. 

6 Transport Measurements 

In this section we present the I{V) vs V curves for the Nb/F/Nb Josephson 
junctions varying the thickness of the F barrier. All these materials show multi- 
ple oscillations of the Josephson critical current with barrier thickness implying 
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Figure 16: Thermal variation of the magnetization for Co, Py and Ni (black 
dots) with the best fitting curves to extrapolate the Curie temperature of Ni 
(571K), Py (800 K) and Co (1200K). Red squares are cooling data for Ni. 

repeated O-tt phase-transitions in the superconducting order parameter. The 
critical current oscillations have been modeled with the clean and dirty limit 
theoretical models, from this analysis we have extrapolated the exchange en- 
ergy and the Fermi velocity of the ferromagnetic barrier. 

6.1 Experimental Data: Critical Current Oscillations 

For several junctions on each chip we have measured the current versus voltage 
from which the critical current (/c) and the normal resistance (Rn) of the 
Josephson junction have been extrapolated. In Fig. [T7| we show an example 
oi I vs V curves for different Co thicknesses, we notice that for each curve we 
have roughly the same resistance, that means the same area of the Josephson 
junctions. This is an indication of the good control on the fabrication of the 
devices with the FIB. 

The IcRn products, as a function of F barrier thickness, exhibit a decaying, 
oscillatory behavior, in agreement with the theoretical predictions. The oscilla- 
tions of IcRn as a function of Ni, Py, Co and Fe thicknesses at 4.2 K are shown 
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Figure 17: I — V characteristics for Nb/Co/Nb Josephson junctions with differ- 
ent Co thickness. 



in Figs, [m (a)-(b) and Figs. [19] (a)-(b) [ffl [HI [H] . 

In the case of Py, the clean limit model, Eq. [5l closely matches the experi- 
mental data up to a thickness of ~ 2 nm, while in the case of Ni the oscillations 
are explained by the clean limit theory up to ~ 7 nm. Above these values a 
better fit is obtained using a formula for a diffusive and high E^^ ferromagnet, 
Eq. [5] For Eq. [5] the fitting parameters are the E^.^,. and the Fermi velocity 
vp. In the case of Eq. [5] the fitting parameters are: vp, the mean free path 
L, the superconducting energy gap A and i^ea;- But L and A are not free pa- 
rameters, because they are fixed by the theoretical predictions. In this way 
the only free parameters are, as for the clean limit equation, the Fermi velocity 
and the exchange energy. Using Eq. [6] we obtain the best fitting values for Ni: 
vpi^i) = 2.8 X 10^ m/s and Lm "^ 7 nm; while for Py Wi^(Py) = 2.2 x 10^ m/s, 
Lpy ~ 2.3 nm and A = 1.3 meV. These values are consistent with the ones 
used in Eq. [5] and elsewhere [32 [52] ; while for the exchange energy we estimate 
i?e2:(Ni) ~ 80 and E^^^y) ~ 201 meV. Ef,x{^\) is close to other reported val- 
ues by photoemission experiments [32] , but smaller than that reported by other 
authors f3Tj. The smaller than expected £'e2,(Ni) is a consequence of impurities 
and possibly interdiffusion of Ni into Nb. Anyway, from the magnetic mea- 
surements the extrapolated value of the Tcuric provides evidence that our Ni is 
of acceptable quality. The i?ea;(Py) is consistent with the expected value and 
is approximately twice that measured in Nb/Py/Nb junctions deposited with 
epitaxial barriers where E^x — 95 meV [22] . 

On the other hand for Co and Fe the thicknesses range all in the clean limit, 
so the experimental data have been modeled with the general formula, Eq. [4] 
and ^1 and ^2 are the two fitting parameters. We can see that the experimental 
data are in good agreement with the theoretical model. 
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Figure 18: Oscillations of the IcRn product as a function of the thickness of the 
Ni (a) and Py (b) barrier. The solid green line is the theoretical fit in agreement 
with Eq. [51 dash-dot blue line is referred to Eq. [HI 

In particular for the Co data, from the theoretical fit shown in Fig. [TOT a) 
we find that the period of oscillations is T = 1.91 nm, hence ^2 ~ 0.30 nm 
and ^1 ~ 3.0 nm. In the clean limit, ^2 — vpfi/'iEex- In this way, known ^2 
and estimating the Fermi velocity of being vp — 2.8 x 10^ m/s, as reported in 
literature, we can calculate the exchange energy of the Co: Eex = Hvp/'^T Ri 
309 meV. 

For Fe, from the theoretical fit we obtain ^1 — 3.8 nm and ^2 — 0.25 nm. 
So the period of the oscillations is T = 1.6 nm. Known ^2 and vp = 1-98 x 10^ 
m/s [33], we can calculate the exchange energy of the Iron: E^x — fivp/2T « 
256 meV. 

To confirm that our oscillations are all in the clean limit (meaning that the 
considered Co and Fe thickness are always smaller than the Co and Fe mean free 
path), we have modeled our data with the simplified formula which holds only in 
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Figure 19: /c-Rat vs tco (a) and tpc (b) with the best fitting theoretical model 
in agreement with Eq. IH dash-dot violet line, and Eq. [5] solid green line. 



this limit, Eq. [5l where in this case Eex and vp are the two fitting parameters. 
From the theoretical fit we obtain Eex{Co) = 309 meV, vf{Co) = 2.8 x 10^ m/s 
and EexiFe) = 256 meV, WF(Fe) = 1.98 x 10^ m/s. We remark that the best 
fits are obtained with exactly the same values as previously reported from Eq. 
[H Both models provide excellent fits to our experimental data showing multiple 
oscillations of the critical current in a tiny (nano-scale) range of thicknesses of 
the Co and Fe barrier. 
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6.2 Estimation of the Mean Free Path 

With a simplified model that is obtained solving the linear Eilenberger equations 
[M] we can estimate the mean free path for Co and Fe. The general formula is: 



tanh 



L 



i 



eff 



ieff io' + L-^ + iCu 



(8) 



where ^e// is the effective decay length given by ^~A = ^f^ + 1^2'^: ^o is the 
Ginzburg-Landau coherence length and ^h is a complex coherence length. In 
the clean limit 1 + LS^^^ > imax{ln(l + L^g"^), ln(L^^^)}. The solution of Eq. 
[S] gives 



?r' = eo'' + i"\^o = 






^^2 — S,H 



vpfi 
2K 



(9) 



and the numerical solution is shown in Fig. 
10t 



]for Co and in Fig. [2T]for Fe. 




2 4 6 8 10 12 14 16 18 20 



L/£. 



Figure 20: Estimation of Co mean free path. The dependence of C2/C1 with 
inverse magnetic length, L/^h, calculated for different ratios of L/^o- Inset: 
inverse decay length, L/^2 = f{L/£,o) for when L/£^h — 0.1 

In the case of Co [16] , following this method we find from Fig. [20] that the 
experimental ratio ^2/^1 — 0.1 corresponds to two inverse magnetic lengths of 
L/£,H — 16.5 and L/£^h — 18.7. By assuming L/S,o — 0.1 and for the estimated 
parameters ^1 ~ 3 nm and ^2 — 0.3 nm we obtained, from the inset in Fig. 
1201 a mean free path L ~ 5 nm. Furthermore to validate the determined mean 
free path of our Co thin film we have estimated ico in a 50 nm thick Co film 
by measuring its resistivity as a function of temperature using the Van der 
Pauw technique. The transport in Co thin films is dominated by free-electron- 
like behavior |35| and hence the maximum mean free path is estimated from 
L = hkp/neC^pi,, where Ue is the electron density for Co and is estimated 
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from the ordinary Hall effect to be 5.8 x 10^* cm 



and pb is the residual 



resistivity. The residual resistance ratio {{pr + Pb)/pb) was measured to be 
~ 1.41 and hence we calculate Lco at 4.2 K to be 10 nm. This verifies that our 
assumption that for Co the oscillations are all in the clean limit is well justified. 
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Figure 21: Estimation of Fe mean free path. ^2/Ci ''■s L/^h, calculated for 
different ratios of L/^o = 0.2. Inset: L/^h vs L/^h to estimate the mean free 
path, L. 

The same method can be done to estimate the mean free path for Fe [T7]. 
In this case the ratio ^2/Ci — 0.06, so supposing L/^o=0.2, we can extrapolate 
from the graphical solution (Fig. [2T|) a value of L/S,h — 29. Considering the 
curve L/£^2 vs L/^h (see inset in Fig. [21]) we obtain a value L/S,2 of about 
25 and we estimate the value of the mean free path of about 6nm. With this 
analysis we can remark that for Co and Fe the condition that all the thicknesses 
are in the clean limit is unambiguously fulfilled. 

6.3 Shapiro Steps and Fraunhofer Pattern 

In Fig. [22] we show an example of typical I vs V characteristics for a sample 
with a barrier of 4.1 nm of Py (a) [21] and 0.8 nm of Fe (b) p7] . 

For these I{V) curves we also show the effect of an applied microwave field 
with an excitation at f — 13.86 GHz and / = 14.01 GHz, respectively. Con- 
stant voltage Shapiro steps appear due to the synchronization of the Josephson 
oscillations on the applied excitation [37]. As expected the steps manifest at 
voltages equal to integer multiples of hf /2e. 

We have looked at the effect of a magnetic field on the maximum supercur- 
rent Imax in our devices. We find that Imax oscillates with applied magnetic 
field, giving rise to a Fraunhofer pattern; however, we also find that Imax, which 
normally corresponds to the central peak of a Fraunhofer pattern, is offset from 
zero applied field to ±Hoffset, which is equal to ±4.9 mT. Fig. [23] shows a 
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Figure 22: A typical I ws V curve of a Josephson junction with 4.1 nm of Py 
(a) and 0.8 nm of Fe (b) (black dash line) and the integer Shapiro steps in the 
voltage-current curve with an excitation at 13.86 GHz (a) and 14.01 GHz (red 
hne). 00 = /i/2e 




Figure 23: A typical Fraunhofer pattern for a device with 2.5 nm thick Fe 
barrier. The central peak is offset by 4.9 mT, which corresponds to a coercitive 
field of the Fe barrier, as confirmed by the hysteresis loop (red line). 



typical Fraunhofer pattern for a Fe barrier device with a barrier thickness of 
~ 2.5 nm. We compared the variation of Imax with applied field to the mag- 
netic hysteresis loop of the same film (measured prior to patterning and device 
fabrication) at 20 K. The offset field is found to correspond approximately to 
the coercive field of the unpatterned film, which is ±Hcoercive — 4.3 mT. The 
central peak is shifted by the coercive field in each direction, which is due to the 
changing magnetization of the ferromagnetic barrier. The side peaks are not 
hysteretic and displaced by the saturation moment of the barrier because the 
hysteresis loop is saturated and the barrier moment is constant for both field 
sweep directions. The coercive field and offset field in the Fraunhofer pattern do 
not exactly agree: firstly, the coercive field is approximated at 20 K; secondly. 
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processing a device in a FIB microscope is likely to harden the Fe magnetic 
domains by virtue of Ga ion implantation. Similar results were measured for 
Co, Ni, and Py. 

6.4 Half Integer Shapiro Steps Close to the O-vr Transition 




f=16.69GHz 



< 




Figure 24: (a) IcRn versus ico compared to the fit of Eq. [5]used as an eye guide 
only; (b) I vs V for ico = 2.1 nm without (dashed blue line) and with (solid 
red line) microwaves at / = 16.69 GHz; this device exhibits a sudden drop in 
differential resistance at V/riipof =0.5. The inset shows the details of the half 
integer Shapiro step. 



In this subsection, we present a study of a single transition from the to 
TT state in Co devices where ico thickness varied from 1.8 to 2.5 nm [15]. The 
Co magnetic dead layer was estimated to be 1.2 nm; although larger than the 
previously estimated dead layer of 0.8 nm for the Co barriers, we found that 
the bulk magnetizations for the two sample sets are similar. Importantly, the 
magnetic data convincingly showed incremental increases in magnetic moment 
with increasing Co thickness. From current-voltage measurements, Ic and Rn 
were extracted so that IcRn could be determined and tracked as a function of 
ico- The characteristic voltage decreases to a small voltage around a mean Co 
barrier thickness of 2.05 nm and then increases, implying a change in phase of 
TT. Each datum point in Fig. I24r a) was measured, and the vertical error bars 
were derived from a combination of estimating Ic and Rj^ from the current- 
voltage curves and from a small noise contribution from the current source. 
Particularly for those devices with the smallest characteristic voltage, there was 
a considerable scatter in the obtained IcRn values. For an eye guide only, we 
have modeled the transition with Eq. [51 assuming Eex{Go) = 309 meV. The 
curve is offset to fit the experimental data. The model ignores any influence of 
a second harmonic term in the current-phase relation. 

In general, the current-phase relation is periodic in 0, the phase difference; 
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however, recent theoretical and experimental works [55] have looked at the pos- 
sibility of observing higher order harmonics, Ig = Id sincj) + Ic2 sin(20), where 
Ic2 ^ Id at the to TT crossover and the second harmonic dominates. When 
Is = Ic2 sin(2(/)), one expects both integer and half-integer Shapiro steps in the 
current-voltage curves at particular microwave frequencies. In our experiment 
[TS], by applying microwaves in the 13-17 GHz range to those devices near the 
transition, we have found that the device with the smallest characteristic volt- 
age and with a Co barrier thickness of 2.1 nm exhibited current steps at both 
half-integer n = 1/2 and integer n — 1 values of V/n(j)of, as plotted in Fig. 
^Ub). These results imply that this device is close to a minimum characteristic 
voltage and provides evidence for a second harmonic in the current-phase rela- 
tion; anyway further investigations must be done before any conclusive remarks 
can be made. 



6.5 Temperature Dependence of the IcRn Product 




V(mV) 

Figure 25: Temperature dependence of the Nb/Co/Nb Josephson junction with 
2.2 nm of Co. 

Further evidence of the oscillatory dependence of the characteristic volt- 
age is given by measurements of its thermal variation. We notice that, for 
each Co thickness, the Josephson junction resistance remains approximately 
constant, slightly increasing when approaching the critical temperature of Nb 
{Tc — 9.2K), as expected; this is shown in Figl25]for tco=2.2 nm. On the other 
hand the critical current, and hence the product IcRn, quickly decreases with 
increasing T and goes to zero at the critical temperature of the Nb. This is 
shown explicitly in Fig. [5S] (a) for different Co thicknesses. 

As for the Co data, in Fig. [57] (a) we show the temperature dependence 
of the IcRn product for different thicknesses of the Fe barrier. We remark 
that also in this case the IcRn product, quickly decreases with increasing the 
temperature [T7] . 
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Figure 26: (a) IcRn vs T for different tco- (b) Osciliations of tfie decay rate 
oi.t{Tq — 4.2K) of Ic'Rn, Eq. (TUl with tlie temperature [stars; the blue line is a 
guide to the eye], as compared with the oscillations of the IcRn product itself 
in Nb/Co/Nb vr junctions. 
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Figure 27: (a) Temperature dependence of the IcRn for different thicknesses of 
the Fe layer, (b) Oscillations of the decay rate aT(To = 4.2K) oi IcRn, Eq.fTOl 
with the temperature [stars; the blue line is a guide to the eye], as compared 
with the oscillations of the IcRn product itself [spheres; the red line is the fit 
given by Eq. 2] in a window of Fe thicknesses. 

It is interesting to notice how the rate at which the critical current drops to 
zero, exhibits a non-monotonic dependence on the Co and Fe layer thickness. We 
can define the relative decay rate ax of the IcRn product with the temperature 
as in Ref. [31], namely 



ariTo) = 



dln[IcRNiT)] 



dT 
dIcRN{T) 



(10) 



T=To 



dT 



T=Ta IcRn (To) ' 



calculated at To = 4.2 K. In Fig. [26] (b) and in Fig. [27] (b) we plot ar as 
a function of ico and ipc, respectively, for consecutive thicknesses; we observe 
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oscillations of ar in phase with the oscillations of the IcRn, in other words the 
smaller IcRn decays to zero with increasing T more slowly than the greater 
IcRn- a satisfactory explanation for this intriguing phenomenon awaits further 
investigation. 

7 Conclusion 

In summary, we have fabricated and investigated Nb/F/Nb Josephson junctions, 
with F standing for a strong ferromagnetic layer of Ni, Py, Co or Fe, varying the 
F barrier thickness. A magnetic dead layer lower than 1.7 nm is extrapolated by 
a linear regression of the thickness dependence of the saturation magnetization. 
From the I{V) vs V curves we have measured the Josephson critical current Ic 
and the normal resistance Rn in order to follow the oscillations of the IcRn 
product as a function of the F thickness. In agreement with the theoretical 
models for the clean and dirty limit, we have fitted the experimental data esti- 
mating the exchange energy and the Fermi velocity of the F barrier (see table 
[2] for a summary of the Eg^ and vp)- Shapiro steps appear at integer multiples 
of the applied voltage. Then, for different thickness of the Co and Fe barrier, 
we have shown that the IcRn product decreases with increasing temperature, 
and in particular the decay rate presents the same oscillatory behavior as the 
critical current. 

With the results of this chapter, we have shown that we are able to produce 
nano-structured Nb/F/Nb 7r-junctions with a high control of the F thickness 
(within an accuracy of 0.2 nm). The low magnetic dead layer gives us the 
possibility to reproduce the transport properties in our heterostructures with 
a small thickness deviation. The estimated E^x is close to that of the bulk F, 
which implies that our films are deposited cleanly with only a small reduction 
in exchange energy. Interfacial roughness and possibly interdiffusion of F into 
Nb is assumed to account for the slightly smaller exchange energy. Moreover, 
in particular for Co and Fe, which exhibit oscillations all in the clean limit, 
the high exchange energy means small period of IcRn oscillations, enabling 
us to obtain the switch from to tt state in a very small range of thicknesses 
(< 3 nm). We can therefore conclude that Co and Fe barrier based Joseph- 
son junctions are viable structures to the development of superconductor-based 
quantum electronic devices. The electrical and magnetic properties of Co and 
Fe are well understood and are routinely used in the magnetics industry: we 
have demonstrated that Nb/Co/Nb and Nb/Fe/Nb hybrids can readily be used 
in controllable two-level quantum information systems. 

In this context, a feasible future development will be to realize pseudo spin 
valve devices using two F layers made of different ferromagnetic materials with 
the different coercive fields separated by a superconducting layer [3S]. In this 
case the spacer layer is relatively thick, and is used to decouple the ferromagnetic 
layers to prevent them from switching at the same field. In such an instance the 
fundamental physics will play with the interlink between these two competing 
orders, and on the practical point of view it could be the starting brick for suit- 
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Table 2: Summary of the exchange energy E^x and the Fermi velocity 
VF extrapolated from theoretical fits of the critical current oscillations 
and the magnetic dead layer to- 



Barrier 


Eex{fne 


V) 


vf{: 


KlO'^ m/s) 


tn(nm) 


Ni 


80 






2.8 


1.7 


Py 


201 






2.8 


0.5 


Co 


309 






2.8 


0.7 


Fe 


256 






1.98 


1.1 



able devices. In this way the magnetization orientation of the F/S/F structure 
can be controlled by a weak magnetic field which by itself is insufficient to de- 
stroy superconductivity, but the small magnetic field enables the softer F layer 
to be aligned with it, while the harder material will not switch. This artificial 
structure allows active control of the magnetic state of the barrier. The mag- 
netoresistance of the pseudo spin valve gives direct access to the information 
about relative orientation of the ferromagnetic layers, and the magnetic state of 
the barrier. The novelty, and strength of such devices, lies in the fact that these 
heterostructures can be realized cither with classical ferromagnetic materials 
and low-Tc superconductors, or with colossal magneto-resistance materials and 
high-Tc cuprates. So what will be the most promising future technology??? 
Will it rely on liquid helium, or nitrogen???? We hope to contribute towards 
providing a proper answer to such dilemmas. On a far-reaching scope, the route 
to room-temperature applications would constitute a more formidable challenge 
for our research. 
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